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Abstract. We report on the identification of 54 embedded clusters around 217 massive protostellar candidates of 
which 34 clusters are new detections. The embedded clusters are identified as stellar surface density enhancements 
in the 2 /rm All Sky Survey (2MASS) data. Because the clusters are all associated with massive stars in their 
earliest evolutionary stage, the clusters should also be in an early stage of evolution. Thus the properties of these 
clusters should reflect properties associated with their formation rather than their evolution. For each cluster, we 
estimate the mass, the morphological type, the photometry and extinction. The clusters in our study, by their 
association with massive protostars and massive outflows, reinstate the notion that massive stars begin to form 
after the first generation of low mass stars have completed their accretion phase. Further, the observed high gas 
densities and accretion rates at the centers of these clusters is consistent with the hypothesis that high mass stars 
form by continuing accretion onto low mass stars. 
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1. Introduction 

Embedded stellar clusters, those clusters that are still sur¬ 
rounded by the molecular clouds in which they formed, are 
the youngest of the stellar clusters. As such, the embed¬ 
ded clusters are of particular interest to understand which 
properties of stellar clusters are related to their origins and 
which are derived from subsequent evolution (Elmegreen 
et al. iMini Lada & Lada 120(1^11 . For example, the mass 
segregation, the concentration of higher mass stars in the 
centers of clusters, that is observed in many optically vis¬ 
ible open clusters is also seen in some of the embedded 
clusters. Because the embedded clusters are too young to 
have undergone significant dynamical evolution, the mass 
segregation must be a property of the process of star for¬ 
mation in the clusters. As another example, open clusters 
exhibit both hierarchical and centrally condensed morpho¬ 
logical types. Observations of both morphological types 
in embedded clusters suggests that the morphology of the 
clusters reflects the morphology of the clouds from which 
the stars formed rather than the dynamical evolution of 
the cluster. Finally, the distribution of stellar masses in 
the embedded clusters ought to be little affected by evo- 
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lution and therefore closest to the initial mass function 
(IMF). 

While the embedded clusters are thought to be among 
the youngest clusters, if we were to identify a class of clus¬ 
ters in which star formation, and therefore the formation 
of the cluster itself, were just beginning, we would po¬ 
tentially be able to address some of the questions as to 
the causes and origins of some of the cluster properties. 
For example, observations of embedded clusters may sug¬ 
gest star formation rather than dynamical evolution as a 
cause of mass segregation, but there remains the question 
of the cause. Is mass segregation a result of the formation 
of massive stars by the collisions and coalescence of lower 
mass stars because collisions will be more common in the 
high stellar density in the center of a cluster (Bonnell et 
al 119981 Testi et al. I1999II ? Is mass segregation a result 
of the formation of massive stars by continuing accretion 
onto existing low mass stars (Beech and Mitalas 119941 
Behrend & Maeder iMin Bernasconi and Maeder ITTMl 
Meynet and Maeder 20001 Keto l2()()'Hll . a process that re¬ 
quires rapid accretion to overcome radiation and thermal 
pressure, and therefore requires dense gas as would be 
found in the center of a dense molecular clouds? Similarly, 
observations of actively forming clusters might potentially 
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address the relationship of the star formation to molecular 
cloud structure. What is the difference in cloud structure 
leading to the hierarchical and centrally condensed mor¬ 
phological types of clusters? Do stars always form first in 
the center of a molecular cloud or in gravitationally col¬ 
lapsing fragments throughout the cloud? Finally, if the 
lower mass stars form first as suggested by several studies 
of open clusters (HerbigESH Stabler IT985II . and required 
by the theories of massive star formation either by coales¬ 
cence or continuing accretion, then observations of clus¬ 
ters in formation may potentially address the origins of 
the IMF. Observations of the stellar mass distribution in 
actively forming clusters may allow the opportunity to see 
mass distributions that are still evolving toward an IMF. 

In this study, we test a hypothesis that a class of ac¬ 
tively forming clusters, the youngest subset of the young 
embedded clusters, may be identified by searching for stel¬ 
lar clusters around previously identified massive protostel- 
lar candidates in isolated molecular clouds. These mas¬ 
sive protostellar candidates, massive stars in their earli¬ 
est stages of formation, have been identified as a class 
of luminous objects having specific IRAS colors that are 
associated with other indicators of massive star forma¬ 
tion such as dense gas and dust, water masers, and ultra¬ 
compact HII (UCHII) regions (Falla et al. 119911 Molinari 
et al. [T!M)i iT!^ i^nnni Sridharan et al. Beuther et 

al l2()()2ajl and associated with isolated molecular clouds. 
These candidates are sources deeply embedded in their 
molecular clouds and therefore probably represent the first 
massive stars to form within the clouds. If there were more 
evolved massive stars in these clouds, then we would ex¬ 
pect their winds, radiation pressure, and supernovae ex¬ 
plosions to have at least partially cleared the region of gas 
and dust revealing perhaps a classic open cluster as would 
typically be found in a more evolved massive star form¬ 
ing region. Thus because the massive protostellar candi¬ 
dates are the first massive stars to form within their host 
molecular clouds, any associated clusters should represent 
clusters in their earliest evolutionary phase. 

Our study is similar in its objectives to a previous sur¬ 
vey of stellar clusters around Herbig Ae/Be stars (Testi 
et al. ITWiniTO . That study also sought to identify very 
young clusters with active star formation using the Herbig 
Ae/Be stars, which are protostars, as indicators of active 
star formation. In our study we use the massive protostel¬ 
lar candidates as indicators of the earliest stages of stellar 
and therefore cluster formation. 

A number of independent near-infrared (NIR) obser¬ 
vations of the regions around several of the massive pro¬ 
tostellar candidates in our target list already show evi¬ 
dence for embedded stellar clusters (references in Table 2). 
In contrast to these previous observations, in this study, 
we undertake a systematic search for embedded clusters 
around all the previously identified high mass protostel¬ 
lar candidates in the lists of Molinari et al. III99()II and 
Sridharan et al. 

We identify the potential clusters as star count density 
enhancements above the mean background level (Lada & 


Lada ll99'^ Carpenter et al. 120001 Ivanov et al. l2002ll using 
the existing K-band observations of the 2MASS database 
(Kleinmann et al. II 99411 . We report on the identification 
of 54 clusters by this technique of which 34 are new de¬ 
tections. We estimate some basic properties that can be 
derived from the J,H,K data in 2MASS. Finally we dis¬ 
cuss some implications of newer theories of high mass star 
formation for the formation and evolution of clusters. 

2. Method 

The 217 targets selected for our study represent the unique 
sources from the combined lists of 163 candidates in 
Molinari et al. dEMI) and 69 in Sridharan et al. 

We searched for clusters within an area of 400" x 400" 
around each of the 217 positions of high mass protostellar 
candidates. 

The 2MASS database contains J,H, and K bands, and 
we used the K-band at 2.2/rm for the initial identification 
of the clusters because K-band suffers the least extinction. 
We constructed stellar surface density maps by spatially 
binning the 2MASS point sources. We selected sources 
with quality flags of A,B,C, and D, excluding lower qual¬ 
ity sources. We used bins of 120"separated by 60"to cover 
the 400"x400"area. Clusters were detected as star count 
density enhancements above the mean background level 
within the sampled area. We used the average noise (cr) in 
the map as a reference and plotted contours of mode-|-2cr 
and above to identify clusters. The mode value represents 
the average star count in the region, and mode-|-2CT re¬ 
veals the enhancement with respect to the average back¬ 
ground. We identified 63 possible clusters for which we 
retrieved data over a larger area of 600" x 600" in all the 
three photometric bands to enable better sampling of the 
background counts and also to construct color-color dia¬ 
grams and color-magnitude diagrams. Based on improved 
sampling of the background that reduced the noise in the 
contour maps, and iterating with different bin sizes, we ex¬ 
cluded some targets that had loose groupings of 5-8 stellar 
sources. After excluding such weak groupings, we identi¬ 
fied 54 clusters all of which contain more than 8 stars, 
although the majority of the clusters show star counts in 
excess of 20-30. We then used a smaller bin-size of 80" (step 
size 40") to reveal relatively small scale structures within 
the detected clusters. The resulting contour maps for each 
of these 54 clusters are presented in Figure. 1. The contour 
levels begin at mode-|-2cr and increase in steps of Icr. 

In order to estimate the number of stars in each cluster 
we must account for the inclusion of background and fore¬ 
ground stars within the defined cluster boundary. To esti¬ 
mate the combined background and foreground stars, we 
sampled for each cluster an adjacent region of I600"^-2500 
"^. The corrected counts are listed in Table. 2 as “true” 
cluster membership. If we assume that within the sensi¬ 
tivity of 2MASS data we might not detect any stars be¬ 
hind the high extinction of a massive star forming cluster, 
then our subtraction of the estimated count of background 
stars would result in an underestimate of the number of 
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true cluster members. Thus our estimate of the numbers of 
stars in a cluster might represent a lower limit to the true 
number in cases where the extinction through the clusters 
is very high. 

2.1. Detection statistics 

We detected 54 clusters out of 217 unique sources which 
implies an overall ~25% detection rate. However, there are 
no clusters detected in the RA range 6hr to 20hr that co¬ 
incides with the Galactic mid-plane although some targets 
were associated with K-band nebulosity. Our technique of 
identifying clusters as stellar density enhancements in the 
2MASS data does not appear suitable to detect clusters 
in the Galactic mid-plane. First, the high extinction in 
the Galactic mid-plane reduces the number of stars that 
we can detect. Second, in the Galactic mid-plane the back¬ 
ground count level is high and star count density enhance¬ 
ments are less significant with respect to the background. 
Excluding the targets in the 6 - 20 hr RA range, the de¬ 
tection rate for clusters around the massive protostellar 
candidates is 60%. 

2.2. Cluster membership and properties 

In Table. 1 we list the radius, distance, and luminosity of 
each cluster. The radius is defined from the area enclosed 
by the mode-|-2(T contours as R = ^arealir. The distance 
and luminosity are generally quoted from their source pa¬ 
pers namely MBGP96 and SriOl. Where available, im¬ 
proved distance estimates are used based on a SIMBAD 
search. In particular, targets beyond the Solar circle are 
relatively less studied making their distance estimates lim¬ 
ited to the values quoted by Wouterloot et al. ifMIHll . A 
histogram of the effective radii of the clusters indicates a 
mean cluster radius of ~lpc. The IRAS peak positions 
and the cluster peak positions coincide in about half of 
the detected sample and do not coincide in the remaining 
half. However, the fraction in which the two peaks coin¬ 
cide, the clusters are noted to be densest and are more 
circularly defined than the rest. 

2.3. Cluster mass estimation 

We estimated the mass of each cluster based on the 
method described by Lada & Lada II2()()3I Hereafter 
LL03). The method incorporates the K-luminosity func¬ 
tion (KLF) models of Muench et al. II2()()2II and the evolu¬ 
tionary tracks of DAntona & Mazitelli fTTOli . The mass 
estimation is done by posing the following question. If 
we were to place the Trapezium cluster at the same dis¬ 
tance as one of the sample clusters and subject it to the 
same dust extinction, how many of its 780 members (see 
Muench et al. 1200211 are we likely to observe? This ques¬ 
tion can be trivially answered if we know the completeness 
limit (limiting magnitude) of our observed data and by 
using DAntona & Mazitelli II1994II tracks along with the 


formula 

Mfe = Kmag - 5 -l- 5 log(distance) - A^. 

We then compute a ratio of the observed number of stars 
to the expected number if our cluster were the Trapezium 
personified. This ratio then merely indicates the ratio of 
the mass of the sample cluster to that of the mass of 
the Trapezium cluster which is assumed to be 413 Mq 
(M uench et al. 1200211 . 

We use the individual extinction values to each cluster 
as derived from the color-color diagrams (see Sec 2.5) of 
each cluster. Further, the 2MASS limiting magnitude in 
K-band of 15mag and an average cluster age of IMyr is 
adopted to compute the masses. The resulting mass es¬ 
timates are listed in Table. 1. Fig. 2 shows the embedded 
cluster mass distribution function (EGMDF) constructed 
using the cluster mass estimates in Table. 1. The solid line 
represents the EGMDF for the sample from our work and 
the dotted line shows the LL03 EGMDF for their sam¬ 
ple of embedded clusters within 2 kpc from the Sun. The 
EGMDF of both samples are similar. 

In Fig. .3 we plot the IRAS luminosity of the targets 
versus the mass of the associated embedded cluster. There 
is a correlation between the mass of the associated em¬ 
bedded cluster and the luminosity of the target massive 
protostar, similar to the correlation found in HAeBe stars 
by Testi et al. H1999II . 

2.4. The morphology of embedded clusters 

In their review of embedded clusters, Lada & Lada (2003) 
suggest that the embedded clusters can be classified into 
two morphological types, hierarchical and centrally con¬ 
densed, and they suggest that these structures may reflect 
the physical processes responsible for cluster formation. 
In particular, the two types may indicate a dominance 
of turbulent (hierarchical) or gravitational (centrally con¬ 
densed) energies. The embedded clusters, and particularly 
the youngest clusters still in formation, offer the opportu¬ 
nity to test this hypothesis against observations. Such a 
test would best be done comparing the structure of the 
stellar distribution of the clusters with the gas distribu¬ 
tion in the clouds. This comparison would require observa¬ 
tions outside the scope of this current research. However, 
we may compare the relative numbers of hierarchical and 
centrally condensed clusters in the two samples of embed¬ 
ded clusters, those identified in our current research as the 
youngest by their association with the massive protostellar 
candidates and those embedded clusters within 2 kpc of 
the Sun, from the list in Lada & Lada (2003). We classify 
the clusters that show single peaks in the stellar density 
maps as centrally symmetric (C-type) and those which 
show more than one peak as hierarchically structured (H- 
type). The results for our sample of clusters (Fig. 1) in¬ 
dicate a ratio H-type/G-type ~0.8. We applied the same 
classification scheme to the 70 embedded clusters of Lada 
& Lada (2003), using stellar density maps produced from 
the 2MASS data by our procedures as described in §2. 
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Fig. 1. Stellar surface density contour plots for the 54 embedded clusters produced with a bin size of 80" (step size 
= 40"). The star symbol represents the position of the luminous IRAS source. Contours begin at mode+2(7 level and 
increase in steps of lu, cr being the average noise measured over the an area of 600"x600"centered on each target. 
The dots represent the point sources with phqual=A,B,C,D catalogued in the 2MASS GATOR database. 
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Fig. 1. contd. Stellar surface density contours of detected clusters 


The nearby embedded clusters have a ratio H-type/C-type 
~ 0.9, essentially the same as their younger counterparts. 

2.5. Photometric analysis 

Once the clusters were detected, we constructed H-K vs 
J-H color-color and H-K vs K color-magnitude (CM) di¬ 
agrams for each of the clusters to investigate the nature 
of the associated point sources. To do this, we first trans¬ 
formed the 2MASS data into the Bessel & Brett system 
using the transformations: 

(H-K)bb = 1.0298 X (H-K)2m.4SS-0.0350 
(J-H)bb = 1.0101 X (J-H) 2MASS- 0.0495 & 

K_b_b = 1.0298 X K 2 MASS-O.O 35 O. 

We plotted all points that fell within the mode-l-lcr con¬ 
tour to allow consideration of even those members that 
fell at the extreme edges of the clusters along with the 
main-sequence dwarf and giant tracks, reddening vectors, 
T-Tauri locus (Meyer et al. I1997II . and the HAeBe locus 
(Lada & Adams I1992II . We dereddened the points that 
fell within the reddening vectors enclosed by the main- 
sequence tracks to the M6-K6 locus in order to find the 
extinction to each point. A histogram of all such extinc¬ 
tion values and an average extinction was constructed for 
each cluster. These values are listed in Table. 2 and are 
used to obtain the corresponding mass estimates for each 
cluster. We also constructed H-K vs K color-magnitude 


diagrams for each cluster and verihed the consistency of 
extinction derived from color-color diagrams. 

2.6. The NIR colors of massive protostars 

The massive protostellar candidates were all initially iden¬ 
tified by their IRAS colors, and the list further refined 
by the association of the candidates with other signs of 
star formation (references in Introduction). We may ask 
whether the NIR colors of any of the clusters or their 
members are distinct. 

For example, in the formation of high mass stars, it has 
been hypothesized that the dense accretion flows that form 
the stars suppresses the formation and later the expan¬ 
sion of an HII region (Walmsley 1995, Keto 2003). Once 
the HII region reaches a size detectable by typical VLA 
observations, the accretion flow will have ended and the 
star should have obtained its main sequence mass. The 
list of massive protostellar candidates contains candidates 
associated with radio continuum emission (HII regions) 
and those without detectable radio continuum. In order to 
identify any differences between the NIR emission of these 
two types of sources, we searched for NIR emission within 
5” of the dust continuum peaks mapped by Molinari et 
al. (2000) and Beuther et al. (2002) using both the 2MASS 
Point Source Catalog (PSC) and Extended Source Catalog 
(XSC). The extended NIR sources are shown as star sym- 
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Fig. 2. Embedded cluster mass distribution function. The 
dotted line represent the ECMDF from LL03 for the em¬ 
bedded clusters within 2 kpc distance. The solid line shows 
the ECMDF for the clusters from this work. 
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luminosity{Log Lsun) 


Fig. 3. A plot of the mass of the associated clusters vs the 
luminosity of the target massive protostar 


bols, and NIR point sources are shown as points in Fig. 4 
plotted along with the main sequence dwarf & giant loci, 
the T-Tauri locus and the HAeBe locus. Reddening vec¬ 
tors assume an interstellar reddening law of R=3.12. The 
sources with no radio continuum (Group A in Fig. 4) dis¬ 
play very red colors, are spread around the region of the 


Group A (Younger) Group B (Older) 



Fig. 4. Color-color diagram of NIR counterparts to the 
1.2 mm peaks of HMPOs. The solid curve represents the 
HAeBe curve by Lada & Adams dMl, thick-solid curve 
represents the main-sequence tracks, the dotted line shows 
the T-Tauri locus from Meyer et al. ITMTIl and the dashed 
lines represent the reddening vectors. The dots represent 
the point sources from GATOR PSC and star symbols 
represent the extended sources from the GATOR XSC. 

main sequence and T-Tauri loci, and extend way beyond 
the limits of the HAeBe locus. In comparison, the older 
sources are limited to the main-sequence region with some 
falling within the HAeBe locus. As far as this test goes, the 
data are consistent with the theory that the massive pro- 
tostellar candidates without associated radio continuum 
are more deeply embedded than those with HII regions. 
The test is limited. The NIR emission, particularly the 
extended emission, is not necessarily associated with any 
massive protostellar candidate. The NIR emission is from 
regions of very high extinction that may affect the J,H, 
and K magnitudes. Spectroscopic observations would be 
more reliable in identifying individual sources. However, 
the J,H, and K data available in the 2MASS survey in 
combination with the continuum emission from dust may 
be sufficient to suggest the locations of the massive pro¬ 
tostellar candidates without suggesting a correspondence 
with individual stars. 

3. Discussion 

3.1. Implications for the relative ages of low and high 
mass stars 

If one accepts the hypothesis that the high mass proto¬ 
stellar candidates identified by the IRAS satellite can be 
localized by the observed combination of the dust contin¬ 
uum and NIR, then there is evidence for mass segregation 
in our sample of the youngest clusters still in formation. 
Furthermore if we accept the hypothesis that the mas¬ 
sive protostellar candidates are the first massive stars to 
form in a cluster, the presence of the cluster itself, which is 
identified primarily through its more numerous lower mass 
members, suggests that the low mass stars form before 
their high mass counterparts. This hypothesis has been 
in the literature for quite some time (e.g. Herbig [T9f)2|l . 
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The clusters in our study, by their association with mas¬ 
sive protostars, reiterate the suggestion with slightly dif¬ 
ferent reasoning than presented earlier (e.g. Hillenbrand 
et al. Ilbbdl Testi et al. I1090II . The difference is that, in 
our sample the sources are relatively more massive than 
typical HAeBe stars, ~40% of the detected clusters are 
also associated with massive outflows (see Table. 2) and 
many sources lack any signihcant HII regions. Therefore, 
the massive stars are still in accreting phases as evidenced 
by their outflows while the low mass members of the as¬ 
sociated clusters have finished accreting and emerged into 
the Class I or II phases as evidenced by their 2^m appear¬ 
ance. Therefore our result places a stronger constraint that 
the massive stars begin to form at least after a few 10^ yrs 
after the hrst low mass stars in these embedded clusters 
were born. 

3.2. Implications for the formation mechanism of 
massive stars 

The previous study of stellar clusters associated with 
Herbig Ae/Be stars by Testi et al. II1997I suggested 

that the stellar spectral type of the Herbig Ae/Be star 
may be correlated with the richness (stellar density) of the 
cluster. This leads to the hypothesis that dynamical inter¬ 
action (collisions) between low mass stars may be required 
to form a high mass star (Testi et al. II 99911 . That hypoth¬ 
esis was motivated in part by the high stellar densities 
observed. Our study does not contradict this suggestion, 
but our results also suggest an alternative hypothesis. 

In our sample, those clusters in which the IRAS/mm 
peaks and the cluster peaks coincide most closely are also 
the densest clusters both in terms of stellar density as 
well as molecular gas density. The very high gas densi¬ 
ties in the centers of these clusters, inferred from obser¬ 
vations of maser and thermal molecular line emission and 
radio continuum emission are consistent with the hypoth¬ 
esis that the high mass stars form by continuing accretion 
onto lower mass stars. The hypothesis of continuing accre¬ 
tion requires very massive accretion flows, perhaps up to 
10“^ Mq yr“^ to form the earliest type stars. (KetolHnSll. 
These massive flows are more likely found in very dense 
gas that our study, in conjunction with previous observa¬ 
tions, locates at the centers of the clusters where the most 
massive stars are forming. 

There are three principle reasons why the the hypoth¬ 
esis of continuing accretion requires a very high rate of 
accretion to form an early type star. First, the rate of ac¬ 
cretion must be high enough to keep a star from evolving 
off the main sequence before it has reached a high mass 
(KetoEISJ- Second, the momentum of the accretion flow 
Mv must be high enough to overcome the high radiation 
pressure of the stellar luminosity fKahn lT974l Wolfire & 
Cassinelli II987I Keto & Wood ITOObll . Third, the density 
of the gas around the star must be high enough to keep 
the HII region around the star small so that the accretion 
flow is not reversed by the thermal pressure of the ionized 


gas. Specifically the radius of ionization equilibrium must 
be within the distance where the escape velocity from the 
star equals the ionized sound speed lKeto l^()2al2(K)2b|l . 

The hypothesis of continuing accretion has specific im¬ 
plications for our understanding of stellar evolution and 
therefore cluster evolution. It is a fundamental assump¬ 
tion of the standard theory of stellar evolution that once 
the mass of a star is fixed, its subsequent evolution is de¬ 
termined. Stars that continue to gain mass by accretion 
obviously do not have a hxed mass, and their evolutionary 
tracks in color-magnitude space or the Hertzprung-Russel 
diagram are different than for non-accreting stars, and fur¬ 
thermore depend on the rate of accretion. For example, 
the stellar structure calculations of Stabler, Shu & Taam 
Ill980all98()bl [TiMll showed that stars of mass greater 
than about 7 Mq begin hydrogen burning while still in 
the accretion phase and thus do not have a protostellar 
phase equivalent to that in lower mass stars. This result 
has been conhrmed in subsequent calculations (Beech and 
Mitalas IM Behrend & Maeder EHni Bernasconi and 
Maeder fi 9961 Meynet and Maeder f2000|l . Second, contin¬ 
uing accretion will supply fresh hydrogen to the star, and 
if the rate is high enough, the fresh hydrogen will prevent 
the star from evolving off the main sequence. Third for 
stars below about 7 Mq, the accretion will provide fresh 
deuterium to the star and prolong the phase of deuterium 
shell burning, altering the pre-main sequence evolution¬ 
ary tracks. We briefly discuss these differences in stellar 
evolution and their consequences for our understanding of 
cluster evolution. 

The accretion rates and time scales necessary for an in¬ 
termediate mass star to evolve to a high mass star can be 
estimated through stellar structure calculations. Figure. 5 
shows evolutionary tracks for stars with continuous accre¬ 
tion at various rates. The evolutionary tracks were cal¬ 
culated for us by Alessandro Chieffi based on the stellar 
structure models described in Chieffi, Staniero & Safaris 
II1995II . The evolutionary tracks are calculated for accre¬ 
tion rates normalized to 6 x 10“® M© yr“^ for a star of 
unity Mq but with different dependencies on the stellar 
mass in each track. While the mode of accretion in high 
mass star formation is not well known, observations in¬ 
dicate that the rates may scale as a power of the stellar 
mass ( Churchwell lI998l Henning l2()()()l Behrend & Maeder 
iMm . We plot accretion rates scaling with powers of the 
stellar mass from 0.5 to 2.0. A rate proportional to the 
square of the mass is the maximum accretion rate attain¬ 
able in purely spherical accretion iBondi [195211 . 

The time scale for the formation of high mass stars is 
quite short. The time required for a star of 4.0 Mq to reach 
the point where it turns off the main sequence toward the 
giant branch is printed on the figure for each accretion 
rate. The main sequence lifetimes for stars of B to late 
O type are roughly between 0.5 and 1.0 Myr. The most 
massive O stars spend such a short time in the accretion 
phase that observed examples of early O stars in the ac¬ 
cretion phase would be expected to be quite rare. Because 
the accretion rate scales as a power of the stellar mass. 
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Fig. 5. An HR-diagram showing evolutionary tracks for 
massive stars that are gaining mass by accretion. The zero 
age main sequence (ZAMS) line is shown annotated by 
stellar mass in solar masses. The four evolutionary tracks 
are for stars accreting with rates (M)that vary as a power 
law of stellar mass (M), with indexes 0.5,1,1.5 & 2 (bot¬ 
tom to top respectively).tzAMS indicates the time spent on 
the ZAMS line for a 4 Mq star before it deviates off the 
line. The accretion rate at the time of deviation, which is 
also the maximum accretion rate for that track is marked 
below each track. 


these precursors to massive stars spend most of their time 
in the lower mass range where the accretion rate is lower 
and the evolution proceeds more slowly. Thus we would 
not expect to see many O type stars in deeply embedded 
clusters that are at the earliest stages of formation. 

Continuing accretion supplies fresh deuterium to the 
stars allowing deuterium burning in the outer shell of the 
star to continue past the time when the deuterium would 
have been exhausted in a non-accreting star. This keeps 
the stellar radius larger and thus the luminosity higher 
than for non-accreting stars. Thus pre-main sequence evo¬ 
lutionary tracks of accreting stars maintain a higher lumi¬ 
nosity as the pre-main sequence star evolves (Fig. 5; see 
also Behrend & Maeder PdOOlll . At around 4 Mq, the lumi¬ 
nosities of both accreting and non-accreting stars are no 
longer dominated by Deuterium shell burning and their 
luminosities are essentially the same. Thus as accreting 
stars above 4 M 0 gain mass they appear to move up the 
main sequence as defined by the standard theory of non¬ 
accreting stars. If all massive stars in a cluster form by 
continuing accretion onto lower mass stars, then there 
will be no massive pre-main sequence stars and no stars 
above 4 M 0 in the HR diagram on high mass pre-main 
sequence tracks. HR diagrams from observations of clus¬ 
ters are consistent with this hypothesis (NGC 2264, NGC 
6530 Stabler 1985; NGG 6611 Hillenbrand et al. 1993). In 
particular, these clusters show large numbers of stars with 
masses less than 4 M 0 above the main sequence, but few 
stars above the main sequence at higher masses. These 


HR diagrams are of course also consistent with the stellar 
evolution of non-accreting stars. In particular, the absence 
of observed pre-main sequence massive stars could be due 
to the brevity of the pre-main sequence phase in massive 
stars. Nevertheless, the determination of the ages of stars 
and clusters by the location of the stars on the HR diagram 
already has some inherent uncertainty (Stabler 1985), and 
the possibility of continuing accretion adds further un¬ 
certainty to this method of age determination. Thus a 
method to date clusters by their association with proto¬ 
stars as in Testi et al. or by their association with 

massive protostellar candidates as in this study should be 
useful in distinguishing between the different hypotheses 
of accreting and non-accreting stellar evolution. 

3.3. Observational evidences for continuing accretion 

Finally, we note that several individual studies of mas¬ 
sive protostars in recent years support the conjecture that 
massive stars are most likely formed by continuing ac¬ 
cretion. First, continuing accretion requires that an in¬ 
termediate mass star such as an A or B star forms by 
accretion and should be surrounded by an envelope mass 
much greater than the protostar itself, so that sufficient 
material is available to feed the intermediate mass star to 
build its mass. Second, the accretion rates should be high 
enough of the order of 10“^ M 0 yr“^. Detailed studies of 
sources such as IRAS20126-I-4104 (Gesaroni et al. I1994|l . 
IRAS20293-f3952 (Beuther et al^^, IRAS07427-2400 
(Kumar et al. 12993611 . G31.41-1-0.31 and G24.78-k0.98 
(Beltran et al. 12995(1 have shown the presence of both 
Keplerian disks and large toroids that contain several hun¬ 
dreds to a few thousand solar masses of material. In the 
case of IRAS07427-2400, the massive star is visible even 
at 2/im and is driving a powerful massive outflow. For in¬ 
stance, in the best studied example of IRAS20126-I-4194, 
while a mass of 7 M 0 is estimated by using the observed 
Keplerain velocities, the luminosity derived from a spec¬ 
tral energy distribution curve (^lO"^ L 0 ) is an order 
of magnitude higher than that expected from a single 
7 M 0 Star (Gesaroni et al. 120051) . Further, the mass ac¬ 
cretion rates derived from the above studies are in the 
range of lO”"* M 0 to as high as 10“^ M 0 (Beltran et 
al. Eoni Eoniii- The observational evidences consistently 
suggest a scenario where intermediate mass objects are 
surrounded by massive envelopes with sufficient accretion 
rates to form a massive star through continuing accretion. 
It should be noted that contrary to this scenario, in low 
mass star formation, the envelope mass is greater than the 
protostar mass only in Glass 9 phase {Menv > M*) and 
during Glass I or H phases Menv < M*. 

4. Summary and Conclusions 

We conducted a systematic search for clustering around 
217 candidate HMPOs chosen from the combined lists 
of MBCP96 and Sri92. We used the 2MASS GATOR 
database and the technique of producing stellar surface 
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density contours to detect clusters. We also searched for 
near-infrared counterparts of the 1.2 mm dust continuum 
peaks associated with all candidate HMPOs. 

1) We find 54 embedded clusters associated with 217 
candidate HMPOs indicating a 25% cluster detection rate. 
All targets lying in the Galactic mid-plane did not show 
any clusters, and we attribute this to the insensitivity of 
the 2MASS data to probe into the large extinctions in the 
Galactic mid-plane region. The detection rate for targets 
away from the mid-plane is 60%. 

2) We estimate the mass of each cluster associated with 
massive protostellar candidates and find that the embed¬ 
ded cluster mass distribution function is similar to that 
found in a sample of all embedded clusters withing 2 kpc 
of the Sun. 

3) Approximately equal numbers of clusters associated 
with massive protostellar candidates are found to have 
hierarchical as centrally condensed structures. 

4) In about half of the detected sample, the cluster 
peaks and the IRAS/mm peaks coincide very well indi¬ 
cating a positive identification of the massive protostar 
with the infrared visible embedded clusters. This fraction 
of clusters also display the highest densities and more cir¬ 
cular morphology in the entire sample. 

5) One hypothesis of the formation of massive stars 
by continuing accretion is that the younger stars will not 
be associated with radio free-free emission from HII re¬ 
gions. We find that the near IR colors of sources near 
the dust continuum peaks and that do not show free-free 
emission are redder than those that are associated with 
ultra-compact HII regions. 

6) The data are consistent with the hypothesis that 
massive stars form by continuing accretion, but the data 
do not discriminate against the hypothesis that massive 
stars form by the collisions of lower mass stars. 
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Table 1. Properties of embedded clusters associated with candidate HMPOs 


EC"- 

No. 

IRAS 

Name 

d L Cluster Members 

kpc Log Lq mode+2(T true 

Effec Radius 
" pc 

Afc : 

mag 

Mstellar Outflow'’ 
Ms 

1* 

00117+6412 

1.8 

3.1 

26 

13 

66 

0.57 

0.3 

18 

Y 

2* 

00420+5530 

4.3 

4.1 

71 

38 

96 

3.6 

0.6 

380 

Y 

3* 

03211+5446 

4.54 

4.5 

131 

41 

101 

2.21 

0.6 

471 

- 

4* 

04034+5116 

3.98 

4.2 

50 

16 

112 

2.16 

0.3 

100 

- 

5 

05137+3919 

11.5 

4.6 

33 

7 

116 

6.08 

0.3 

1132 

Y 

6* 

05168+3634 

6.08 

4.4 

48 

18 

115 

3.37 

0.4 

355 

Y 

7 

05274+3345 

1.55 

3.6 

48 

29 

103 

0.77 

0.5 

38 

Y 

8 

05345+3157 

2.6 

3.6 

95 

36 

126 

1.09 

0.6 

122 

Y 

9 

05358+3543 

1.8 

3.8 

53 

30 

130 

1.13 

0.5 

48 

Y 

10* 

05373+2349 

2.4 

3.3 

11 

7 

65 

0.88 

0.5 

18 

Y 

11* 

05382+3547 

16.5 

5.7 

98 

60 

152 

18.76 

0.8 

- 

- 

12 

05490+2658 

2.1 

3.5 

95 

43 

122 

1.23 

0.6 

94 

Y 

13 

05553+1631 

2.5 

3.8 

80 

22 

104 

1.24 

0.4 

56 

Y 

14 

06056+2131 

1.5 

3.8 

132 

74 

132 

0.95 

0.8 

111 

Y 

15 

06061+2151 

2.0 

1.4 

105 

49 

145 

0.07 

1.0 

144 

- 

16* 

06063+2040 

4.52 

4.9 

202 

114 

121 

2.64 

0.7 

1474 

- 

17 

06068+2030 

1.5 

4.7 

68 

34 

96 

0.69 

0.6 

45 

- 

18 

06103+1523 

4.3 

4.63 

50 

24 

118 

2.64 

0.3 

211 

- 

19 

06104+1524 

4.7 

4.68 

50 

24 

124 

2.81 

0.3 

217 

- 

20* 

06105+1756 

4.5 

3.38 

70 

35 

114 

1.87 

1. 

314 

- 

21 

06155+2319 

1.6 

1.3 

127 

71 

132 

1.02 

0.6 

102 

- 

22 

06291+0421 

1.96 

1.96 

38 

20 

113 

1.07 

0.5 

36 

- 

23 

06308+0402 

2.02 

3.9 

47 

17 

96 

0.94 

0.6 

34 

Y 

24* 

06382+0939 

0.76 

2.2 

72 

38 

125 

0.46 

0.5 

28 

- 

25 

06584-0852 

4.48 

3.9 

32 

16 

108 

2.33 

0.4 

145 

Y 

26 

20028+2903 

1.55 

3.7 

37 

10 

72 

1.29 

0.7 

15 

- 

27 

20050+2720 

0.73 

2.6 

57 

28 

105 

0.37 

0.7 

21 

Y 

28* 

20081+2720 

0.7 

2.5 

36 

12 

78 

0.26 

1.0 

9 

- 

29* 

20205+3948 

4.5 

4.5 

44 

17 

74 

1.61 

0.6 

190 

- 

30 

20286+4105 

3.72 

3.6 

44 

7 

58 

1.00 

0.8 

63 

Y 

31* 

20319+3958 

1.6 

3.8 

32 

9 

68 

0.53 

1. 

17 

- 

32* 

20406+4555 

11.9 

5.1 

64 

5 

28 

0.69 

0.5 

1766 

- 

33* 

20444+4629 

2.42 

3.5 

81 

17 

57 

0.96 

0.3 

37 

N 

34* 

21046+5110 

0.59 

2.4 

62 

24 

85 

0.99 

0.6 

17 

- 

35* 

21078+5211 

1.49 

4.1 

113 

50 

113 

0.81 

0.6 

66 

Y 

36* 

21202+5157 

6.78 

4.8 

44 

3 

84 

2.77 

0.5 

98 

- 

37* 

21526+5728 

8.11 

4.4 

24 

4 

53 

1.12 

0.4 

202 

- 

38* 

21548+5747 

7.1 

4.4 

15 

5 

38 

1.3 

0.3 

140 

- 

39 

22134+5834 

2.6 

4.1 

52 

26 

104 

1.31 

0.4 

72 

Y 

40* 

22147+5948 

7.29 

4.4 

59 

5 

106 

2.26 

0.4 

178 

- 

41* 

22172+5549 

2.87 

3.7 

171 

35 

126 

1.74 

0.3 

109 

Y 

42* 

22272+6358 

1.23 

3.3 

23 

11 

74 

0.44 

0.4 

11 

- 

43* 

22305+5803 

5.4 

4.1 

36 

3 

83 

2.16 

0.4 

43 

Y 

44* 

22506+5944 

5.7 

4.3 

38 

15 

91 

2.52 

0.5 

280 

Y 

45* 

22551+6221 

0.7 

3.2 

229 

129 

152 

0.51 

0.8 

98 

- 

46* 

22570+5912 

5.1 

4.7 

96 

25 

82 

2.01 

0.3 

270 

Y 

47* 

23026+5948 

5.76 

4.2 

59 

22 

86 

2.39 

0.3 

331 

N 

48 

23033+5951 

3.5 

4. 

10 

3 

74 

1.26 

0.3 

14 

Y 

49* 

23146+5954 

4.43 

4.4 

56 

19 

119 

2.55 

0.3 

153 

- 

50* 

23151+5912 

5.7 

5. 

27 

2 

97 

2.69 

0.4 

33 

Y 

51* 

23314+6033 

2.78 

4. 

176 

61 

126 

1.69 

0.6 

239 

Y 

52* 

23448+6010 

2.02 

3.4 

64 

38 

89 

0.86 

0.5 

66 

- 

53* 

23507+6230 

4.28 

4.1 

23 

8 

88 

1.81 

0.4 

66 

- 

54* 

23545+6508 

0.8 

3.6 

92 

39 

145 

0.56 

0.7 

31 

- 


* denotes newly detected embedded clusters. The remaining clusters were known previously as listed below 

(5)Ishii et al. I2()()2j i7.12.17.21.23') Carpenter et al. 119941 (8)Chen et al. I1999j (18,26,30) Dutra & Bica 120(111 (9) Porras et 

al l2000l (13,14,15) Carpenter et al. 119951 (22) Phelps & Lada fTMTI : (25) Ivanov et al. 120021 (39,48) Kumar et al. 120021 


Outflow identifications from Beuther et al l20fl2b1 Zhang et al l2005l and Wu et al. 120041 



















